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•!. Introduction 



Magnetically-ordered quasicrystals: Enumeration of 
spin groups and calculation of magnetic selection rules 

Ron Lifshitz* and Shahar Even-Dar Mandel 

School of Physics and Astronomy Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv Uni- 
versity, Tel Aviv 69978, Israel. Correspondence e-mail: ronlif@post.tau.ac.il 

We provide the details of the theory of magnetic symmetry in quasicrystals, which 
has previously only been outlined. We develop a practical formalism for the enu- 
meration of spin point groups and spin space groups, and for the calculation of 
selection rules for neutron scattering experiments. We demonstrate the formalism 
using the simple, yet non-trivial, example of magnetically-ordered octagonal qua- 
sicrystals in two dimensions. In a companion paper [Acta Crystallographica AXX 
(2003) xxx-xxx], we provide complete results for octagonal quasicrystals in three 
dimensions. 
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yhe discovery of Zn-Mg-RE icosahedral quasicrystals, where 
RE stands for one of the rare-earth metals Y, Gd, Tb, Dy, Ho, 
and Er (Niikura et al, 1994; Tsai et al, 1994), and initial indi- 
cations that some of these quasicrystals may possess long-range 
I ]antiferromagnetic order (Charrier et al, 1997; Charrier et al., 
1998), have generated increased interest in the nature of mag- 
netic order in quasicrystals (Fukamichi, 1999, and references 
therein). Subsequent measurements (Islam et al., 1998; Sato 
et al., 1998), followed by ongoing vigorous research (Fisher 
et al, 1999; Sato et al, 1999; Fisher et al, 2000; Sato 
'et al, 2000a; Sato et al, 2000^; Dolinsek et al, 2001; Kramer 
et al, 2002), have shown that only short-range spin correla- 
tions exist in these particular quasicrystals, and therefore the 
existence of long-range magnetic order in real quasicrystals re- 
mains an unresolved question of great interest. The recent dis- 
covery of cadmium-based quasicrystals (Guo et al, 2000fl; Guo 
et al, lOQQb; Tsai et al, 2000) and initial studies of mag- 
netic order in the Cd-Mg-Tb icosahedral quasicrystal (Sato 
et al, 2002), may provide some insight into this question. The- 
oretical models that deal with magnetism on quasicrystals — 
purely geometrical models (Niizeki, 1990a; Niizeki, \99Qb; 
Lifshitz, 1995; Lifshitz, 2000) as well as physical ones, such 
as the Ising model (Grimm & Baake, 1997; Matsuo et al, 
2000; Matsuo et al, 2002), the XY model (Hermisson, 2000), 
the Heisenberg model (Wessel et al, 2003), and the Hubbard 
model (Jagannathan & Schulz, 1997; Hida, 2001) — are known 
to exhibit long-range magnetic order. This is despite initial in- 
tuition that aperiodicity necessarily induces geometrical frustra- 
tion and is therefore inconsistent with having magnetic order 

It is therefore clear — ^both from a theoretical and an experi- 
mental standpoint — that there is a need for a theoretical classi- 
fication of all types of quasiperiodic long-range magnetic order 
that are allowed by symmetry. Such a symmetry-based classi- 
fication of quasiperiodic magnetic structures, combined with a 
calculation of the selection rules, imposed by magnetic sym- 

' We usually consider 3-dimensional magnetic moments, or spins, in a rf-dimensional crystal, where = 2 or 3, and therefore take S(r) to be a 3-component field. 
If necessary, one can generalize to spins of arbitrary dimension. 



metry, will be of great assistance in the study of real mag- 
netic quasicrystals, if and when they are discovered. It will also 
offer valuable guidance in the search for these novel materi- 
als (Lifshitz, 2000). In this paper we provide the details of the 
extension to quasicrystals (Lifshitz, 1998) of Litvin and Ope- 
chowski's theory of spin groups (Litvin, 1973; Litvin & Ope- 
chowski, 1974; Litvin, 1977). In section 2. we explain how 
spin groups are used to describe the magnetic symmetry of pe- 
riodic as well as quasiperiodic crystals. In section 3. we give 
the group theoretic details of the classification of the different 
groups into appropriate equivalence classes. In section 4. we de- 
velop the formalism required for the actual enumeration of spin 
groups, and in section 5. we derive the relations between the 
magnetic symmetry of a crystal and the selection rules it im- 
poses on neutron diffraction experiments. In section 6. we treat 
the case of octagonal symmetry in two dimensions as a peda- 
gogical example. In a companion paper (Even-Dar Mandel & 
Lifshitz, 2003) we perform the actual detailed enumeration of 
three-dimensional octagonal spin groups which turns out to be 
surprisingly rich. Complete enumeration for the other common 
quasicrystals without explicit details will follow in future pub- 
lications. 

2. Using spin groups to describe the symmetry of 
magnetically-ordered crystals 

A (i-dimensional magnetically-ordered crystal, whether peri- 
odic or aperiodic, is most directly described by its spin den- 
sity field S(r). This field is a 3-component real- valued function, 
transforming like an axial vector under 0(3) and changing sign 
under time inversion. One may think of this function as defin- 
ing a set of classical magnetic moments, or spins, on the atomic 
sites of the material.^ For quasiperiodic crystals the spin den- 
sity field may be expressed as a Fourier sum with a countable 
infinity of wave vectors 



S(r) = ^S(k) 



,ik r 



(1) 
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The set L of all integral linear combinations of the wave vectors 
in (1) is called the magnetic lattice, and is characterized among 
other things by a rank D, and by a lattice point group Gc- Its 
rank D is the smallest number of wave vectors needed to gener- 
ate it by integral linear combinations. For quasiperiodic crystals, 
by definition, the rank is finite. For the special case of periodic 
crystals the rank is equal to the dimension d of physical space. 
The set of (proper or improper) rotations, that when applied to 
the origin of Fourier space leave the magnetic lattice invariant, 
is the lattice point group Gl, also called the holohedry. 

The theory of magnetic symmetry in quasiperiodic crys- 
tals (Lifshitz, 1998), is a reformulation of Litvin and Ope- 
chowski's theory of spin space groups (Litvin, 1973; Litvin 
& Opechowski, 1974; Litvin, 1977). Their theory, which is 
applicable to periodic crystals, is extended to quasiperiodic 
crystals by following the ideas of Rokhsar, Wright, and Mer- 
min's "Fourier-space approach" to crystallography (Rokhsar 
et at., 1988a; Rokhsar et ai, 1988^).^ At the heart of this ap- 
proach is a redefinition of the concept of point-group symme- 
try which enables one to treat quasicrystals directly in physi- 
cal space, as opposed to the alternative "superspace approach" 
(Janssen et at., 1992). The key to this redefinition is the ob- 
servation that point-group rotations (proper or improper), when 
applied to a quasiperiodic crystal, do not leave the crystal in- 
variant but rather take it into one that contains the same spatial 
distributions of bounded structures of arbitrary size. 

This generalized notion of symmetry, termed indistinguisha- 
bility, is captured by requiring that any symmetry operation of 
the magnetic crystal leave invariant all spatially-averaged auto- 
correlation functions of its spin density field S(r), for any order 
n and for any choice of components a, G {x, y, z}, 

Cila„(ri,...,r„)= lim i / dr5„,(ri - r) • • •5„„(r„ - r). 

V^oo V Jy 

(2) 

It has been shown (Lifshitz, 1997, in the Appendix) that 
an equivalent statement for the indistinguishability of any two 
quasiperiodic spin density fields, S(r) and S'(r), is that their 
Fourier coefficients are related by 

S'(k) =e2-'x(k)s(k)^ (3) 

where x. called a gauge function, is a real-valued scalar func- 
tion which is linear (modulo integers) on the magnetic lattice L. 
This simply means that 

Vki,k2eL: x(ki+k2)=x(ki) + x(k2), (4) 

where "=" denotes equality modulo integers. 

With this in mind, we define the point group G of a d- 
dimensional magnetic crystal to be the set of operations g in 
0(d) that leave it indistinguishable to within rotations 7 in spin 
space, possibly combined with time inversion.^ Accordingly, 
for every pair {g, 7) there exists a gauge function, (k), called 
a phase function, which satisfies 

Sigk) = e^'^'*?('')7S(k). (5) 

^ For a review, see Mermin (1992fo) or Mermin (1999), for an elementary introc 
' Note that since S(r) is an axial vector field we can restrict 7 to be a proper rot 



In general, as we shall see later, there may be many spin-space 
operations 7 that, when combined with a given physical-space 
rotation g, satisfy the point-group condition (5). We denote 
physical-space rotations by Latin letters and spin-space oper- 
ations by Greek letters. We use a primed Greek letter to explic- 
itly denote the fact that a spin-space rotation is followed by time 
inversion. Thus, the identity rotation in physical-space is e, the 
identity rotation in spin-space is e, and time inversion is e'. Also 
note that we use the same symbol 7 to denote an abstract spin- 
space operation, and to denote the 3 x3 matrix, operating on the 
field S, representing this operation. 

If [g, 7) and (h, rf) both satisfy the point group condition (5), 
then it follows from the equahty 

S([^A]k) = S{g[h\i]), (6) 

that so does [gh, 777). This establishes that the set F of all trans- 
formations 7 forms a group, and the set Gs of all pairs [g, 7), 
satisfying the point-group condition (5), also forms a group. The 
latter is a subgroup of G x F, called the spin point group. We 
shall consider here only finite groups G and F, although in gen- 
eral this need not be the case. The equality (6) further implies 
that the corresponding phase functions, one for each pair in Gs, 
must satisfy the group compatibility condition, 

V(^, 7), % v)eGs: ^];!{k)=^]{hk) + a.^(k). (7) 

Note that successive application of the group compatibihty con- 
dition (7) reveales a relatively simple relation between the phase 
functions of two conjugate elements {g, 7) and {hgh~^ ,r]jr]~^) 
of Gs, 

y{g,i),ih,v)€Gs: ^^;;_";'(/jk)=3.T'(k)+^,^(gk-k). (8) 

A spin space-group, describing the symmetry of a magnetic 
crystal, whether periodic or aperiodic, is thus given by a mag- 
netic lattice L, a spin point group Gs, and a set of phase func- 
tions 'I'g (k), satisfying the group compatibility condition (7). 
We continue to call this a spin space-group even though its 
physical-space part is no longer a subgroup of the Euclidean 
group E{d). Nevertheless, the spin space-group may be given 
an algebraic structure of a group of ordered triplets (g, 7, <I>J) 
in a manner similar to the one shown originally by Rabson et al. 
(1988), and more recently again by Drager & Mermin (1996), in 
the context of ordinary space groups for nonmagnetic crystals. 

In the case of periodic crystals one can show (Mermin, 
1992/7, Eq. (2.18)) that any gauge function 27rx(k), relating two 
indistinguishable spin density fields as in Eq. (3), is necessarily 
of the form k ■ t for some constant translation vector t indepen- 
dent of k, so that S'(r) = S(r + t) and indistinguishability re- 
duces to identity to within a translation. One can then combine 
rotations in physical space and in spin space with translations 

m, see Lifshitz (1996fo). 

1 without any loss of generality. 
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to recover the traditional spin space-groups of periodic crystals, 
containing operations that satisfy 

S(^r)=7S(r + t2'), (9) 

leaving the spin density field identical to what it was. In the 
quasiperiodic case one must retain the general form of OJ(k) 
which is defined only on the magnetic lattice and cannot be hn- 
early extended to arbitrary k. 

3. Classification of spin groups 

The common symmetry properties of different magnetic struc- 
tures become clear only after they are classified into properly 
chosen equivalence classes. We are concerned here with the 
classification of magnetic crystals into Bravais classes (sec- 
tion 3.1.), spin geometric crystal classes (section 3.2.), spin 
arithmetic crystal classes (section 3.3.), and spin space-group 
types (section 3.4.). 

3.1. Bravais classes 

Magnetic crystals, as well as nonmagnetic crystals, are clas- 
sified into Bravais classes according to their lattices of wave 
vectors. Intuitively, two magnetic lattices are in the same Bra- 
vais class if they have the same rank D and point group Gl (to 
within a spatial reorientation) and if one can 'interpolate' be- 
tween them with a sequence of lattices, all with the same point 
group and rank. Stated more formally, as presented by Drager 
& Mermin (1996), we say that two magnetic lattices L and L' 
belong to the same Bravais class if: 

1. The two lattices are isomorphic as abelian groups, i.e. 
there is a one to one mapping, denoted by a prime ('), 
from L onto L' 




satisfying 

(ki + kz)' = k'l + k^ ; (11) 

2. The corresponding lattice point groups Gl and G'l are 
conjugate subgroups of 0{d), 

G'l = tGlt-' , (12) 

for some proper li-dimensional rotation r; and 

3. The isomorphism (10) between the lattices preserves the 
actions of their point groups, namely 

(^k)' = ^'k' , (13) 

where ^ = rgr^^. 

Since the classification of magnetic lattices for magnetic crys- 
tals is the same as the classification of ordinary lattices for 
nomagnetic crystals we shall not expand on this issue further 
but rather refer the interested reader to previous discussions on 
the matter (Rokhsar et al., 1987; Mermin et al., 1987; Mermin 
et al, 1990; Mermin, 1992/;; Mermin & Lifshitz, 1992; Mer- 
min, 1992a; Lifshitz, 1996^; Drager & Mermin, 1996). 



3.2. Spin geometric crystal classes 

When we say that two magnetic crystals "have the same spin 
point group" we normally mean that they belong to the same 
equivalence class of spin point groups, called a spin geometric 
crystal class. We say that two spin point groups Gs and G'g are 
in the same spin geometric crystal class if they are conjugate 
subgroups of 0{d) x [S0{3) x 1'], where 1' is the time inversion 
group, containing the identity e and the time inversion operation 
e'. This simply means that 

G's=ir,a)Gs{r,a)-\ (14) 

for some physical-space rotation r G 0(d), and some spin- 
space operation a E 5(9(3) x 1'. The effect of these rotations 
on the spin point group Gs is to reorient its symmetry axes both 
in physical space and in spin space. 

3.3. Spin arithmetic crystal classes 

The concept of a spin arithmetic crystal class is used to dis- 
tinguish between magnetic crystals which have equivalent mag- 
netic lattices and equivalent spin point groups but differ in the 
manner in which the lattice and the spin point group are com- 
bined. Two magnetic crystals belong to the same spin arithmetic 
crystal class if their magnetic lattices are in the same Bravais 
class, their spin point groups are in the same spin geometric 
crystal class, and it is possible to choose the lattice isomorphism 
(10) such that the proper rotation r used in (12) to establish the 
lattice equivalence is the same rotation used in (14) to establish 
the spin point group equivalence. 

3.4. Spin space-group types 

The finer classification of crystals in a given spin arithmetic 
crystal class into spin space-group types is an organization of 
sets of phase functions into equivalence classes according to 
two criteria: 

1 . Two indistinguishable magnetic crystals S and S', related 
as in (3) by a gauge function x, should clearly belong to 
the same spin space-group type. Such crystals are neces- 
sarily in the same spin arithmetic crystal class but the sets 
of phase functions <I> and <I>' used to describe their space 
groups may, in general, be different. It follows directly 
from (3) and from the point group condition (5) that two 
such sets of phase functions are related by 

07(k)=a.J(k) + x(gk-k), (15) 

for every (g, 7) in the spin point group and every k in 
the magnetic lattice. We call two sets of phase functions 
that describe indistinguishable spin density fields gauge- 
equivalent and equation (15), converting <I> into <!>', a 
gauge transformation. The freedom to choose a gauge x 
by which to transform the Fourier coefficients S(k) of the 
spin density field and all the phase functions <I>, describ- 
ing a given magnetic crystal, is associated in the case of 
periodic magnetic crystals with the freedom one has in 
choosing the real-space origin about which all the point 
group operations are applied. 
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2. Two distinguishable magnetic crystals S and S', whose 
spin space-groups are given by magnetic lattices L and 
L', spin point groups Gs and G^, and sets of phase func- 
tions <t> and <!>', have the same spin space-group type if 
they are in the same spin arithmetic crystal class, and if, 
to within a gauge transformation (15), the lattice isomor- 
phism (10) taking every k e L into a k' e L' preserves 
the values of all the phase functions 

4>7;(k')^<i>j(k), (16) 

where g' = rgr~^ and 7' = a^a~^ . Two sets of phase 
functions that are related in this way are called scale- 
equivalent. This nomenclature reflects the fact that the 
lattice isomorphism (10) used to relate the two magnetic 
lattices may often be achieved by rescaling the wave vec- 
tors of one lattice into those of the other. 

4. Enumeration of spin groups 

The task of enumerating spin groups is limited to the enumer- 
ation of the distinct types of spin point groups and spin space- 
groups. This is because the classification of magnetic lattices 
into Bravais classes, as well as the determination of all distinct 
relative orientations of point groups G with respect to these 
lattices, giving rise to different arithmetic crystal classes, are 
the same as for nonmagnetic crystals, and therefore need not 
be enumerated again. The enumeration of possible spin point 
groups and spin space-groups is greatly simpUfied if one first 
lists all the general conatraints these groups must obey due to 
their algebraic structure. We list below the general constraints 
on the spin point group Gs (section 4.1.), discuss the conse- 
quences of these constraints on the group of spin-space opera- 
tions r (section 4.2.), describe a particularly interesting connec- 
tion between a certain subgroup of T and the magnetic lattice L 
(section 4.3.), and then outline the sequence of steps taken in 
the enumeration of spin groups (section 4.4.). 

4.1. Structure of the spin point group Gs 

The algebraic structure of the spin point group Gs is severely 
constrained by the point group condition (5) as described by the 
five statements below. Proofs for the first four statements can be 
found in the review on color symmetry (Lifshitz, 1997, Section 
IV.A.) as they apply equally to the structure of the color point 
group of a colored crystal. 

1 . The set of real-space operations associated with the spin- 
space identity e forms a normal subgroup of G, called Ge. 

Note that as a special case of Eq. (8) the phase functions 
of conjugate elements of G^ are related by 

e G„ {h, 7))&Gs: 3.^,,-: (/tk)=a.|(k)+a.^(^k-k). 

(17) 

2. The set of spin-space operations paired with the real- 
space identity e forms a normal subgroup of F, called 

the lattice spin group F^. In the special case of peri- 
odic crystals, the elements of F^ are spin-space operations 



that, when combined with translations, leave the mag- 
netic crystal invariant. 

Again, as a special case of Eq. (8) the phase functions of 
conjugate elements of F^ are related by 

V7 e F„ {h, ri)eGs: 3.,^^""' (/jk)=<J.7(k). (18) 

3. The lattice spin group F^. is abelian. 

4. The pairs in Gs associate all the elements of each coset of 
Ge with all the elements of a single corresponding coset 
of Fg. This correspondence between cosets is an isomor- 
phism between the quotient groups G/Ge and F/F,.. 

5. If two phase functions 'I>7'(k) and <I>p(k), associated 
with the lattice spin group F^, are identical on all wave 
vectors then 71 =72- 

Proof : From the point group condition (5) we obtain 

VkGL: 7iS(k) = 72S(k), (19) 

implying that the two operations have the same effect on 

all the spin density fields whose symmetry is described 
by this paricular spin space group type, and are therefore 
identical. 

4.2. Consequences for F and 

The lattice spin group F,, is severely constrained by being an 

abelian subgroup of S0{3) x 1'. Namely, it can have no more 
than a single axis of n-fold symmetry with n > 2. This implies 
that the possible lattice spin groups F^ are the ones listed in the 
first column of Table 1 . 

The fact that the lattice spin group F^ is a normal subgroup 
of F implies that F cannot contain any rotation a G 50(3) for 
which o-FgO-"' 7^ Fg. One can easily verify that the possible su- 
pergroups F for each lattice spin group F^ are the ones listed in 
the second column of Table 1. 

4.3. Relation between the magnetic lattice L and the lattice 
spin group F^. 

We have already mentioned that in the special case of peri- 
odic crystals, the lattice spin group F^, is the set of all spin-space 
operations that, when combined with real-space translations, 
leave the magnetic crystal invariant. It should be of no surprise 
then that in the quasiperiodic case there should remain an inti- 
mate relation between the lattice spin group F^, and the magnetic 
lattice L. We describe this relation here without proof, which 
can be found in the review on color symmetry (Lifshitz, 1997, 
Section IV.C.) where a similar relation exists between the lattice 
L and lattice color group of a colored crystal. 

Recalling that the lattice L is itself an abelian group under 
the addition of wave vectors, one can show that it necessarily 
contains a sublattice Lq, invariant under the point group G, for 
which the quotient group L/Lq is isomorphic to the lattice spin 
group Fg. This isomorphism is established through the prop- 
erties of the phase functions OJ(k) associated with all ele- 
ments 7 of the lattice spin group. In particular, the sublattice 
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Lq is defined as the set of wave vectors k for which the phases 
<J>J(k)=0, for all elements 7 of the lattice spin group. Further- 
more, the relation (18) between phase functions of conjugate 
elements of Eg ensures that the isomorphism between L/Lq and 
Eg is invariant under all elements {h, t]) of the spin point group. 
In other words, if the isomorphism maps a particular wave vec- 
tor k to a particular spin operation 7, then for every {h, rj) in Gs, 
the wave vector hk is mapped to rijri~^. 

This relation between the lattice spin group and the magnetic 
lattice not only imposes a severe constraint on the possible lat- 
tice spin groups but also provides an additional method to cal- 
culate the phase functions (k). One of two alternative ap- 
proaches can be taken to enumerate the allowed combinations 
of Fe and F: 

1. For each type of lattice spin group F^, listed in Table 1, 
see whether there exists an invariant sublattice Lq of 
L giving a modular lattice L/Lq isomorphic to F^, and 
whether the possible extensions of F^. into supergroups F, 
also listed in Table 1, allow the isomorphism to be invari- 
ant under the spin point group. 

2. For each type of lattice spin group F^. and its possible 
extensions into supergroups F, listed in Table 1, simply 
try to solve all the group compatibilty conditions (7) im- 
posed on the phase functions 'i>]{k), associated with the 
elements of F^ and the wave vectors of L. If a solution 
exists then F^ is a possible lattice spin group, otherwise it 
is not. 

It should be emphasized that, either way, the possible combina- 
tions of F,, and F, and therefore the possible types of spin point 
groups, cannot be determined independently of the choice of 
magnetic lattice L. 

4.4. Enumeration steps 

The enumeration of spin point groups and spin space-groups 
consists of a sequence of steps which are listed schematically in 
the flowchart of Fig. 1. We shall illustrate the whole process in 
section 6. by enumerating, as an example, all the 2-dimensional 
octagonal spin point groups and spin space-groups. 

One begins by choosing a lattice L from any of the known 
Bravais classes. One then picks any point group G, compati- 
ble with L, and lists all its normal subgroups along with the 
corresponding quotient groups G/G^. One then chooses one of 
the normal subgroups and calculates, using one of the two 
approaches described in the previous section, all allowed com- 
binations of F and F^. such that the quotient group F/F^ is iso- 
morphic to G/Gg. One then pairs the cosets of G^ in G with the 
cosets of F,, in F in all distinct ways. After checking for equiv- 
alence, as described in section 3.2., one arrives at a list of the 
distinct types of spin point groups. 

For each spin point group one then looks for all solutions to 
the group compatibility conditions (7) not already considered 
above. These solutions are organized into gauge-equivalence 
and scale-equivalence classes, as described in section 3.4., 
yielding the distinct spin space-group types. Because phase 
functions are linear on the lattice L [Eq. (4)] it is sufficient 



to specify their values on a chosen set of D wave vectors that 
primitively generate the lattice. Also, it is sufficient to specify 
the phase functions only for a small set of operations (g, 7) that 
generate the spin point group. All other phase functions can be 
determined through the group compatibility condition. Further- 
more, one can greatly simpUfy the calculation of phase func- 
tions by making a judicious choice of gauge prior to solving 
the group compatibilty conditions, rather than solving the group 
compatibihty conditions and only then organizing the solutions 
into gauge-equivalence classes. 

5. Calculation of magnetic selection rules 

Magnetic selections rules, or symmetry-imposed constraints on 
the form of the spin density field, offer one of the most direct 
experimental observations of the detailed magnetic symmetry 
of a magnetic crystal. In elastic neutron scattering experiments, 
every wave vector k in L is a candidate for a magnetic Bragg 
peak, whose intensity is given by (Izyumov & Ozerov, 1970) 

/(k)oc|S(k)|2-|k.S(k)|2, (20) 

where k is the scattering wave vector and k is a unit vector in 
its direction. It has been shown (Lifshitz, 1996a) that, under 
generic circumstances, there can be only three reasons for not 
observing a magnetic Bragg peak at k even though k is in L\ 
(a) The intensity /(k) ^ but is too weak to be detected in the 
actual experiment; (b) The intensity /(k) = because S(k) is 
parallel to k; and (c) The intensity /(k) = because magnetic 
selection rules require the Fourier coefficient S(k) to vanish. 
Selection rules that lead to a full extinction of a Bragg peak are 
the most dramatic and easiest to observe experimentally. Other 
types of selection rules {e.g. that lead to an extinction of one of 
the components of S(k), or to a nontrivial relation between the 
components of S(k)) are harder to observe. 

We calculate the symmetry-imposed constraints on S(k), for 
any given wave vector k e L, by examining all spin point-group 
operations {g, 7) for which gk = k. These elements form a sub- 
group of the spin point group which we call the little spin group 
of k, G5. For elements {g, 7) of G^, the point-group condition 
(5) can be rewritten as 

7S(k) = e-2'^*?«S(k). (21) 

This implies that every Fourier coefficient S (k) is required to be 
a simultaneous eigenvector of all spin-space operations 7 in the 
little spin group of k, with the eigenvalues given by the corre- 
sponding phase functions. If a non-trivial 3-dimensional vector 
satisfying Eq. (21) does not exist then S(k) will necessarily van- 
ish. It should be noted that the phase values in Eq. (21) are inde- 
pendent of the choice of gauge (15), and are therefore uniquely 
determined by the spin space-group type of the crystal. 

The process of determining the form of the simultaneous 
eigenvector S(k) is greatly simplified if one makes the fol- 
lowing observation. Due to the group compatibility condition 
(7) the set of eigenvalues in Eq. (21) for all the elements 
{g, 7) G G^ forms a 1 -dimensional representation of that group. 
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Spin space-group symmetry thus requires the Fourier coeffi- 
cient S(k) to transform under a particular 1 -dimensional rep- 
resentation of the spin-space operations in the little spin group 
of k. We also independently know that S (k) transforms under 
spin-space rotations as a 3-dimensional axial vector, changing 
its sign under time inversion. It is therefore enough to check 
whether the particular 1 -dimensional representation, dictated 
by the spin space-group, is contained within the 3-dimensional 
axial- vector representation. If it is not, then S(k) must vanish; 
if it is, then S(k) must lie in the subspace of spin space trans- 
forming under this 1 -dimensional representation. 

6. Octagonal spin groups in two dimensions — An 
example 

To demonstrate the ideas presented in this paper, we enumer- 
ate the octagonal spin groups in two dimensions, and calculate 
the magnetic selection rules that arise for each spin space-group 
type. We choose to treat the octagonal crystal system because it 
is the most interesting example for a magnetic quasicrystal in 
two-dimensions. The reason for this is twofold: First of all, as 
for nonmagnetic 2-dimensional crystals, only when the order of 
synnmetry is a power of 2 is it possible to have space groups 
with nonsymmorphic operations; Secondly, only when the or- 
der of symmetry is a power of 2 is it possible to have simple 
imtiferromagnetic order (Niizeki, 1990a; Niizeki, 1990^?; Lif- 
shitz, 1997; Lifshitz, 2000). 

Only partial enumerations of spin groups on quasicrystals ex- 
ist to date. Decagonal spin point groups and spin space-group 
types in two dimensions have been listed by Lifshitz (1995) 
without providing much detail regarding the enumeration pro- 
cess. All possible lattice spin groups for icosahedral qua- 
sicrystals have been tabulated (Lifshitz, 1998) along with the se- 
lection rules that they impose, but a complete enumeration of all 
icosahedral spin groups was not given. This is therefore, the first 
complete and rigorous enumeration of spin groups and selec- 
tion rules for a quasiperiodic crystal system in any dimension. 
In a companion paper (Even-Dar Mandel & Lifshitz, 2003) we 
enumerate the octagonal spin groups in three dimensions, and 
in future publications we intend to treat all the other common 
quasiperiodic crystal systems, though we shall probably not in- 
clude the full details of the calculation. 

6.1. Two-dimensional octagonal point groups and Bravais 
classes 

The lowest rank D that a two-dimensional octagonal lat- 
tice can have is 4. There is just a single Bravais class of two- 
dimensional rank-4 octagonal lattices (Mermin et ai, 1987). All 
lattices in this two-dimensional Bravais class contain an 8-fold 
star of wave vectors of equal length, separated by angles of j 
(as shown in Fig. 2), of which four, labeled b'^^ . . .b'''^ can 
be taken as integrally-independent lattice generating vectors. 
The lattice point group Gl is 8mm, generated by an eightfold 
rotation rg and either a mirror of type m, which contains one 
of the generating vectors and its negative, or a mirror of type 
m', which lies between two of the generating vectors. The two- 
dimensional point groups G to be considered in the enumeration 



are 8mm and its subgroup 8. There is only a single way to orient 
the two point groups with respect to the lattice, so there is just 
a single spin arithmetic crystal class for each spin geometric 
crystal class. 

6.2. Enumeration of spin point groups 

We begin by listing in the first columns of Table 2 all normal 
subgroups of the point groups G = Smm and 8, and the re- 
sulting quotient groups G/G^- Note that the two subgroups 2mm 
and m of the point group 8mm are not normal and therefore do 
not appear in the Table. 

As generators of the spin point groups we take the generators 
of G (rs and m for G = 8mm, and rg alone for G = 8), and 
combine each one with a representative spin-space operation 
from the coset of with which it is paired. We denote the spin- 
space operation paired with by S and the operation paired 
with m by n. When (or m) are in G^ we take S (or /i) to be 
e. The constraints on the operations S and /i, due to the isomor- 
phism between G/G^ and F/F^,, are summarized in the fourth 
column of Table 2. To the generators (rg, (5) and (m, /i) we add 
as many generators of the form (e, 7,) as required, where 7, are 
the generators of F^ (three at the most). Although this set of 
spin point-group generators may, in general, be overcomplete, 
it is the most convenient set to take. 

6.3. Calculation of possible F and 

We use the group compatibility conditions (7) on the phase 
functions ^>J(k), associated with elements in the lattice spin 
group Fg, in order to calculate the possible combinations of F 
and Fg. 

We first note, from inspection of Table 2, that no quotient 
group G/Gc contains an operation of order 3. This implies, 
among other things, that F/F^ cannot contain such an operation 
and therefore the extensions of the orthorhombic lattice spin 
groups Fe, listed in the third row of Table 1, into supergroups 
F cannot be cubic — they can be tetragonal at most. This then 
implies that for any possible combination of F and F^, 

V7 (E F„ 5 e F : (5^7(5-2 = 7. (22) 

With this relation at hand we can proceed to prove the fol- 
lowing short lemmas: 

1. The lattice spin group F^ contains no more than 3 ele- 
ments 7 7^ e, all of which are of order 2. 

Proof : Let (5 G F be the operation paired with rg in the 
spin point group. The relation (22) together with Eq. ( 1 8), 
relating phase functions of conjugate elements in F^, 
yield 

07 (b(') )=3>f T^"' (r^b(') ) (r^b(') ) . (23) 
Thus, for any 7 G F^ 

(24) 

and 

4>J(-b«)=4>7(b(')) => 4>7(b('))=0 or ^. (25) 
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The last result (due to the linearity of the phase func- 
tion), impUes through the group compatibiUty condition 
that 4>7 (b^'^)=0, and therefore, that 7^ = e, or that 7 is 
an operation of order 2. It also implies that each of the 
phases a and /3 in (24) can be either or 1 /2, but they 
cannot both be if 7 7^ e. Thus, there can be no more 
than 3 operations in other than the identity. 

2. Only a single element 7 ^ e in the lattice spin group Eg 
commutes with the operation ^ e F, paired with r% in the 
spin point group. 

Proof : If 7 7^ e in commutes with 5, then the rela- 
tion (18), between phase functions of conjugate elements 
of 7e, implies that 

a>7(b('))=OfT*"' (r8b('))=<J>7(r8b(')). (26) 

Thus, 7 is necessarily the operation whose phase function 
is given by (24) with a=/3=l/2. 

These lemmas, together with the facts that G/G^ can be no 
bigger than a group isomorphic to 8mm, and that the order of 
the operation 6 paired with is no bigger than 8 (proven in 
section 6.4.1. below), narrow down the possible combinations 
of F and F^., listed in Table I, to the ones hsted in Table 3. 

6.4. Enumeration of spin space-group types 

We now turn to the enumeration of spin space-group types 
by calculating the possible values of the phase functions for the 
generators (rg, 6) and (m, fi) on the star of generating vectors 

6.4.1. The phase function for {n, 5). As in the case of reg- 
ular space groups for nonmagnetic crystals (Rokhsar et ai, 
1988fc; Rabson etal., 1991), there is a gauge in which the phase 
function <Pf^ (k) = on the whole lattice. This can be shown 
by starting with arbitrary values for the phase function and 
performing a gauge transformation (15) with the gauge function 

X(b») = i< |^|]b(^)j , i = l...4, (27) 

where b<^'' = -b^^""^' for j = 5,6,7,8. The change to <^f^ 
caused by this gauge transformation exactly cancels it 

A<(b(')) = X (r8b(') - b(')) =^<(b('+4) - b(')) 

= -<(b»), (28) 

so that after the gauge transformation <&f^ (k) = for all wave 
vectors k. Note that this implies, through the group compatibil- 
ity condition (7), that (k)=0, so that 6^ = e, imposing an 
additional restriction on the group F, as indicated in Table 3. 

6.4.2. The phase function for (m, /u). When G = 8mm we 
need to calculate the additional phase function 4>^(k), associ- 
ated with the second point-group generator (m, /i). The gen- 
erating relations that contribute to the determination of this 



phase function are (m, /i)^ = (e, /i^) and (rg, S){m, n){r$, S) = 
{m,5fj,5). Applying the group compatibility condition (7) to 
these relations, in the gauge where (k)=0, yields 



a)^'(b('))=3)^(mb(''-Fb«), 
«I>f,^''(b«)=<I>^(r8bW). 



(29) 



(30) 



We shall first determine the value of the phase 3>^(b^'') us- 
ing Eq. (29), and then use Eq. (30) to infer the values of 
on the remaining three generating vectors. We start by noting 
that iJ,^ is an operation in F^ which is the square of an opera- 
tion in F. Inspection of all the possibilities, listed in Table 3, 
reveals that only two operations, 2^ and e, satisfy this condi- 
tion. Furthermore, if m is the mirror that leaves b^'^ invariant, 
then application of Eq. (29) to b^^' which is perpendicular to m 
(OTb(3) = -b(3)) yields 



«Df (b(^')=«I)^(mb^'' + b^'')=0. 



(3) J- h(3)^ 



(31) 



This implies that /U^ cannot be 2^, because 4>^= has the value 5 
on all lattice generating vectors. Therefore, /x^ must be equal to 
e. Application of Eq. (29) to b^'^ now yields 



0^2<I>^(b('))^4>^(b('))=0ori 



(32) 



and application of Eq. (29) to b'^^ and b'"*^ shows that 
<I>^(b'^^)=<I>^(b but provides no further information re- 
garding the actual values of these phases. 

Next, we examine Eq. (30), which can be rephrased (using 
the group compatibility condition (7)) as 



3)^(bW) + "'^''^(bW)=^^(r8bW). 



(33) 



The value of <t>,^, on b''' determines the values of <I>^ on the 
remaining generating vectors through some phase function, as- 
sociated with an element of F^. Note that fjr^SfiS is an opera- 
tion in Fe which is the product of two operations, /i^'(5/i and 
S, that are conjugate in F. Inspection of Table 3 shows that if 
the product of any two conjugate operation in F is in F^., then 
this product is necessarily either 2^ or the identity e. Substituting 
the values (b^'^ ) = 0000 and <^l' (b^'^ ) = i i i i, we conclude 
that 

^^:(b«)JT"Vr = (34) 

""^ ' [OiOiorioiO if5ii6 = ii2i. 

Thus, there are two spin space groups for each 2-dimensional 
octagonal spin point group with G = 8mOT. 

6.5. Spin group tables 

The resulting 2-dimensional octagonal spin point groups and 
spin space groups are listed in Table 4 for G = 8, and in Table 
5 for G = 8mm, using the following notation: 

Each line in the Tables represents one or more spin point 
groups and their associated spin space groups. The spin point 
groups are given by their generators, listed in the fifth column 
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of each Table. The first four columns provide the group theo- 
retic structure of the spin point group by listing the lattice spin 
group Fe, the normal subgroup Ge of G, paired with Eg, the quo- 
tient group G/Ge, and the full group of spin-space rotations F, 
satisfying the requirement that G/G^ — r/Fg. We use stars and 
daggers to denote optional primes on elements of T (i.e. the ap- 
plication of time inversion after a spin space rotation). If two 
operations in F can be independently primed or unprimed we 
use a star for the first and a dagger for the second. For example, 
the symbol 2*2*2 stands for the two possible groups F = 222 
and 2'2'2, whereas the symbol 4*2^2*^ stands for four distinct 
groups, F = 422, 4'22', 42'2', and 4'2'2. 

To list the spin space groups for each spin point group we 
must specify the values of the phase functions for all the spin 
point-group generators on the four generating vectors of the lat- 
tice. The phase functions <t>J for generators of the form (e, 7) 
are already listed in Table 3 and are not repeated in Tables 4 and 
5. The phase function <I>f^ is zero everywhere due to the choice 
of gauge, and is therefore also not listed in Tables 4 and 5. The 
two possible values of the phase function <t>,^, when the point 
group is 8mm, which according to Eq. (34) depend on the value 
of 5iJ,5, are Usted in the sixth column of Table 5.lf S^S — 
we write "0; i" to indicate that a>^(bW)=0000 or When 
SfiS ^ /i we write "A" to indicate that 4>,^,(b^'')=050i or jO^Q, 
alternating its value from one generating vector to the next. 

In the last column of each Table we give a unique sym- 
bol for each spin space group, based on the familiar Interna- 
tional (Hermann-Mauguin) symbols for the regular (nonmag- 
netic) space groups. To incorporate all the spin space-group 
information we augment the regular symbol in the following 
ways: (1) The symbol for the lattice spin group F^ is added in 
parentheses immediately after the regular space group symbol, 
unless Fe = 1 . (2) In the case of 2-dimensional octagonal spin 
space groups, the values of the phase functions associated with 
the elements of F<, are unique and therefore need not be listed. 
In general, one can encode these phase functions by indicat- 
ing the sublattice Lq (for which L/Lq is isomorphic to F^, as 
described in section 4.3.) as a subscript of the magnetic lattice 
symbol P. (3) To each generator of the point group G we add a 
superscript listing an operation from the coset of F,, with which 
it is paired (if that operation is e we omit it, if it is e' we simply 
add a prime, we use stars and daggers, as described above, to 
denote multiple groups, and we omit the axis about which ro- 
tations are performed if it is the f-axis). (4) The value of phase 
function <J>,^,, when the point group is 8mm, is encoded by re- 
placing the secondary m by a ^ (as in the International symbols) 
when ^^(b^'') = 5555, and by adding a subscript a (for "al- 
ternating") so that nia indicates that <J>^(b^'^) = O^Oj and ba 
indicates that a)^(b(')) = ^0^0. 

6.6. Selection rules due to 

Operations in G which impose selection rules for neutron 
diffraction are those that leave some lattice vectors invariant. In 
two dimensions these can be the identity e, which leaves all lat- 
tice vectors invariant or mirror lines that leave all vectors along 
them invariant. We first consider the selection rules that arise 



from operations (e, 7) where 7 e Fg, and therefore apply to all 
lattice vectors, expressed in terms of the four generating vectors 
as k = mb^^) -I- n2b(^) + nsb*^' + uM'^l 

6.6.1. Selection rules for Fe = 2, 2', 1'. Denoting the gener- 
ator of Fg by 7, the phases <i>] (b'^'^)=^j^j in all of these cases. 
This implies through the eigenvalue relation (21) that the form 
of S(k) depends on 7 and on the parity of ^ n, as follows 

7S(k) =e-™S"'S(k). (35) 

Namely, whenever ^ n,- is even the phase in Eq. (35) vanishes 
and S (k) must be invariant under the operation 7; and whenever 
^ n, is odd the phase is in and S(k) must change its sign un- 
der 7. The consequences for the three possible operations 7 are 
summarized in Table 6. 

6.6.2. Selection rules for Te = 222, 22' 2'. Here F^. is gen- 
erated by (e,2|) and (e,2p, with phase functions given by 

(bW)=OiO^ and^f (b(')) = iOiO. The eigenvalue relations 
(21) for the two generators are 

2*S{k) = e-''^(''^+«)S(k), (36) 
2|S(k) = e-''^('"+"')S(k), (37) 

so that S(k) remains invariant (changes its sign) under 2| if 
n2 + nn is even (odd); and remains invariant (changes its sign) 
under 2| if ni + ^3 is even (odd). These results are summarized 
in Table 7 for the two possible T^s. 

6.7. Selection rules on mirror lines 

In addition to the selection rules arising from F<, there are also 
selection rules that occur when k lies on one of the mirror lines 
and is therefore invariant under reflection through that particular 
mirror. In this case the eigenvalue equation (21) imposes further 
restrictions on the Fourier coefficients of the spin density field. 

Vectors lying along the mirror m„ that leaves the generating 
vector b^'' invariant, have the general form 

k,- = n,bW + /,(b<'-'^ -I- b''+')), i = 1,2, 3,4, (38) 

where all indices are taken modulo 8, and b'-'' = — b'-' "*' for 
j = 5, 6, 7, 8. Selection rules along mi, which is the mirror m 
used to generate the point group (see Figure 2), are determined 
by the equation 

MS(ki) = e-2''^"'*"('''"'S(ki), (39) 

where we have used flie fact [Eq. (34)] that ^^(b*^') - 
<I>^(b''*')=0. Therefore, the form of S(ki) depends on /X, on 
the parity of rii, and on the phase «I>^(b^^^) as follows: If ni is 
odd and <I>,^,(b^'') = 5 then S(ki) must change its sign under 
11, otherwise S(ki) must remain invariant under /x. 

To obtain the selection rules for vectors lying along the re- 
maining mirrors m,- (/ = 2, 3,4), we use successive applications 
of the symmetry operation (r, 6) to the result (39) for mi. Since 
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mki = ki, it follows from relation (8), between phase functions 
of conjugate operations, that 

6'-'iid'-'S{ki) = e-2'™'*i;(''<")s(k,.). (40) 

Thus, in general, for a vector k,, given by Eq. (38) and lying 
along the mirror m;, the form of S(k,) must satisfy the fol- 
lowing requirement: If n,- is odd and «I>^(b'^') = | then S(k,) 
must change its sign under otherwise S(k,) must re- 

main invariant under (5'^'/i(5'^'. Note that <I>,^,(b''') = for 
spin space groups of type P%^m^^m{Te) and P^^m'^miTe)^ and 
that <I>^(b^^^) = 5 for spin space groups of type P8^b^m{re) 
and PS^b^m{re). Also note that in most cases F is abelian, so 
5'~^IJ,d^~' = jjL, and the selection rules take a much simpler 
form. These results are sunmiarized in the second colunm of 
Table 8. 

Vectors lying along the mirror m-, which is between the gen- 
erating vector b^'^ and b'''''^\ have the general form 

k; =n,(b«-fb('+'))+/,(b('-''+b('+2)), i = 1,2,3,4, (41) 

again with all indices taken modulo 8, and b^^' = — b'^""*^ 
for i = 5,6,7,8. Using the group compatibility condition for 
the relation m[ = r^nii (see Figure 2), in the gauge where 
3>f^(k) =0, yields 

<^(k'i) ^ O^(k'i) ^ (a.^(bW) +Oi;(b(2))) , (42) 

where we have used the fact [Eq. (34)] that 3)^(b(^')=4>^(b(')) 
and (-b*"*^ ) =4>^ (b<^' ) . Therefore, the selection rules for m[ 
are deternoined by the equation 

5/xS(k;) = e-2-(".+'.)K(b'")+<i.j;(bp'))s(k'i), (43) 

and again, by successive rotations (r, 8), we obtain the selection 
rules for the remaining mirrors m[, 

<5'M(5i-'s(k;) = e-2-(«.+'.)(*i:(b'")+*-(b'^'))s(k/)_ (44) 

Thus, in general, for a vector k-, given by Eq. (41) and lying 
along the mirror m[, the form of S(k;) must satisfy the follow- 
ing requirement: If n,- + U is odd and 4>^(b*^') + 0^(b*^') = \ 
then S(k-) must change its sign under d'iJ,d^~', otherwise S(k-) 
must remain invariant under 5'n5^~\ Note that <l>^(b'^') + 
«I>,^,(b'^^-') = for spin space groups of type P%^m'^m{re) and 
P2,^b>'m{Te), and that 4>i5(b<'') + *i;(b'^') = i for spin space 
groups of type P%^m'^m{re) and P?,^b'^m{re)- Also note that in 
most cases Y is abelian, so = Sfi, and the selection rules 

take a much simpler form. These results are summarized in the 
third column of Table 8. 
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Table 1 

Possible lattice spin groups Ve and their extensions into the full groups F of spin-space operations. All possible F^'s are listed in the first column. The second column 
shows the constraints on F imposed by the fact that is a normal subgroup of F. The integer k is arbitrary. 





F 




F C 50(3) xl' 


n; n';nl' 


F C (te)221' 


222;2'2'2' 


F C 4321' 


2j2'2' 


F C 4221' 



Table 2 

Normal subgroups Ge of the point groups G = Smm and 8. The resulting quotient group G/Gc is represented in the third column by a point group, isomorphic to it. 
Constraints on the spin-space orperations S and fi, paired with the generators rg and m of G are listed in the fourth column. In each line the first power of 5 that is in 
F(. is given, fi^ is always in F^, therefore we only note whether £ F,,. If 5 or /i are in F^ they are chosen as e. 



G Ge 


G/Ge 


Constraints 


8mm Smm 


1 


5 = E, n = e 


8 


m 


5 = £,H fVe 


4mm 


2 


52 GFe,// = e 


4m'm' 


2 


5 = Al ^ Fe 


4 


2mm 




2 


4mm 




1 


Smm 




8 8 


1 


5 = e 


4 


2 


5^ GF, 


2 


4 


5" GF, 


1 


8 


5« GF, 
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Tables 

Possible lattice spin groups Te and their extensions into the full groups T of spin-space operations, compatible with the 2-diniensional rank-4 octagonal lattice. All 
possible r^'s are Usted in the first column and constraints on the possible supergroups F are listed in the second column. The phase functions for the generators of 
Ve are Usted in the third column. 





Constraints on F 


1 


F^ C F C 8221' 


1' 


Fe C F C 8221' 


2 


Fe C F C 8221' 


2' 


F<. C F C 8221' 


222 


422 C F C 4221' 


2j2'2' 


42'2' C F C 4221' 



Phase functions for generators of F^ 



N/A 



(bW) : 



1111 

2 2 2 2 

ilii 

, 2 2 2 2 

\" ' — 2122 

2iVhWA =oioi; 

2w._f,u _ ^Inl . 



(b«) 
<I>^*(b«) 



(b«) = inl 



iOiO 



Table 4 

Two-dimensional octagonal spin point groups and spin space group types with point group G = 8. The phase functions is zero everywhere by choice of gauge. 
The values of the phase functions for 7 £ on the lattice generating vectors are Usted in Table 3. The symbols for the spin space groups are Usted in the 
rightmost column using the notation described in section 6.5.. 







G/G, 


F 


Generators 


Space Groups 


1 


8 


1 


1 




PS 




4 


2 


2* 


{'•8,21) 


P82* 








1' 




P8' 




2 


4 


4* 


in, 41) 


PS'** 




1 


8 


8* 


{'•8,81) 


PS8* 










('■8,8|*) 


PS«'* 


2 


8 


1 


2 


('■8,e), («,2f) 


PS(2) 




4 


2 


4* 


{r8,4|),(e,2,-) 


PS'** (2) 








2*2*2 


('-8:2*)(«,2j) 


PS^I (2) 








21' 


('-8,e'){«,2j) 


PS' (2) 




2 


4 


8* 


(r8,8|){e,2,-) 


P88* (2) 


2' 


8 


1 


2' 


('•8,^){<^,20 


PS(2') 




4 


2 


2'2'2 


(r8,2|)(e,20 


ps^; (2') 








21' 


{'•8,e'){^,20 


PS' (2') 




2 


4 


41' 


('-8:4i)(f,29 


PS-* (2') 




1 


8 


81' 


('-8,8?)(<^,20 


PS* (2') 


1' 


8 


1 


1' 


{r8,e){e, e') 


PS(l') 




4 


2 


21' 


('-8.2j)(«, e') 


PS2(1') 




2 


4 


41' 


('-8,4j)(«,e') 


PS'*(1') 




1 


8 


81' 


{r8,8,-){e,e') 


PS«(1') 










{rg.,&j)(m,e') 


ps«'(r) 


222 


4 


2 


4*22* 




PS**' (222) 
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Table 5: Two-dimensional octagonal spin point groups and spin space 
group types with point group G = %mm. The phase functions <t>^j^ is zero 
everywhere by choice of gauge. The values of the phase functions <t>J 
for 7 e Eg on the lattice generating vectors are listed in Table 3. The 
possible values of the phase function <I>^ are listed in the sixth column 
using the notation described in section 6.5.. The symbols for the spin 
space groups are listed in the rightmost column using the notation which 
is also described in section 6.5.. 





Gc 


G/G, 


1 


Generators 




Spin Space Group Types 


1 




1 


1 







1 


P%mm\ 


P'dhiii 




8 


2 


2* 


(r8,e)(m,2j;) 





\ 


PEm^ m; 


PSb^'m 








1 


(r8,e)(m, e ) 





\ 

1. 


PSm'm; 


PSb'm 




4-mm 


2 


2* 


(r8,2j)(m,e) 





1 


f 8^ mm; 


Pi^* bm 








1 


(r8, e j(m, e) 





1 

?. 


PS' mm; 


PS'bm 




All 


2 


2* 


(r8,2|)(m,2|) 





1 

?. 


r>o2* 2* 

P8^ m"^ m; 


rtrt')* 1 9* 








1 


(r8,e )(m, e ) 





1 

?, 


PS m m; 


P%'h'm 




4 


2mm 


2 2'2 ' 


(r8,2|)(m,2i) 





1 

2 












ml 
21 


(r8,2j(w,2^-) 





1 
?, 


Pq2 ^2 * 


r^o t? til 















1 
2 














(rc, f'^fm ?*\ 
l'^8i t Jl_m, Zi j 





1 

7, 


PQ/™2* ™. 


/^o £7 171 








zL*9t9*t 
T- z, z. 







1 

?. 


pc4* ^2'^ J-.. 


PC4* l2^ 




1 
i 


Off Iff I 


8*7t7*t 
o z, z. 







1 

?. 


po8* ™2l9. 


po8* /,2l 















1 

?. 


PS^"m'^2; 


P8^"b^^m 


2 


8WOT 


1 


2 


(r8,e),(m,e)(e,2^-) 





1 

?. 


P8mm(2y, 


PSbm{2) 




8 


m 


2*2*2 


(r8,e)(m,2*)(e,2.) 





1 

?. 


P8m2?m(2); 


PSb^^m{2) 








21' 


(r8,e)(m, e')(e,2f) 





1 

?. 


P8m'm(2); 


PWm{2) 




4mm 


2 


2*2*2 


(r8,2|)(m,e)(e,2j) 





1 
?, 


f82'-mm(2); 


P8^^bm{2) 








21' 


(r8,e')(m,e)(e,2j) 





1 

2 


PS'mm{2y, 


PS'bm{2) 








4* 


(r8,4*)(m,e)(e,2,) 


A 


P%'^' mam{2); 


P8^'bam{2) 




4m'm' 


2 


2*2*2 


(r8,2*)(m,2|)(e,2,-) 




Pi^:m'^im{2); 


PS^''b^^m{2) 








21' 


(r8,e')(m,e')(e,2j) 




P8'm'm(2); 


P8'b'm{2) 




4 


2mm 


4*2t2*t 


(r8,4|)(m,2|)(^;,2,) 




P8'**m2*m(2); 


PS^'0m{2) 










(r8,2|)(m,2|,)(e,2,-) 


A 


P82ima^m(2); 


P?Abfm{2) 








2221' 


(r8,e')(m,2*)(e,2.) 


o;| 


P8'm2*m(2); 


P8'b^'m{2) 










(r8,2|)(m,e')(e,2,:) 




P82fm'm(2); 


P8^^b'm{2) 










(r8,2*)(m,2|')(e,2,) 




P82Jm2"m(2); 


P8^' b^'^' m(2) 




2 


4mm 


8*2t2*t 


(r8,8|)(m,2l)(e,2,-) 


0;1 


P8^*m2'-m(2); 


PS^'0mi2) 



continued on next page 
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Table 5: continued 





Ge 


G/Ge 


i 


Generators 


m 


Spin Space Group Types 


2' 




1 


2' 


(r8,e)(OT,e)(e,29 





J 

2 


P8mm(2'); 


PSbm{2') 




8 


m 


2'22' 


(r8,e)(OT,2jE)(e,20 





\ 
2 


P8/n^'m(2'); 


PSb^'m{2') 








21 


(r8,e)(m,e')(e,2i) 





\ 

J 


P8m'm(2'); 


PWm{2') 




4mm 


2 


2'22' 


{rs,2i){m,e)[e,2^} 





2 


P8?-'mm{2'); 


PS'^'bmi2') 








21 


{rn,f'){m,f){e.2',) 





I 

2 


F8'mm(2'); 


P%'bm{2') 




All 

4m m 


2 


2'22' 


[i-^.T,)[m.l^)[e.l'.) 





1 


P?,in-' m-' ni(2'): 


P?,in-- h-- m(2') 








21 


[n,e'){m,f'){e,2,) 





? 


P8'm'm(2'); 


P8'b'm{2') 




4 


2mm 


2221' 


(r8,e')(m,2|)(e,2;) 





2 


P8 m^'m{2'); 


P8 b^'m[2') 










(r8,2^-)(m,e )(e,2,) 





I 


P2/-'m'm{2'); 


Pi^'b'm{2') 




2 


4mm 


4221' 


(r8,4|)(OT,2^-)(e,29 





1 

?. 


P8^ m'^'m{2'); 


PS m 'm{2 ) 




1 


8 mm 


8221 


(r8,8j)(OT,2jE)(e,2j) 





1 

?. 


P8 m 'm{2 j; 


Pj!, b 'myl ) 


i 


dimm 


i 


i 


(r8,e)(m, e)(e, e') 





1 

?. 


P8mm(l'); 


P%bm{V) 




8 


m 


21' 


(''8,e)(OT,2f)(e,e') 





1 

2 


P8m2m(l'); 


P8b^m{V) 




4mm 


2 


21' 


('"8,2z)(w,e)(e,e') 





1 

2 


f82mm(l'); 


P82fom(l') 




4m'm' 


2 


21' 


(r8,2,-)(m,2j)(e,e') 





1 
2 


PSWmCl'); 


P82fc2^(l') 




4 


2mm 


2221' 


(r8,22)(m,2.f)(e, e') 





1 

2 


P82m2^m(l'); 


P8V*m(l') 




2 


4mm 


4221' 


(r8,4.)(m,2^.)(e,e') 





1 
2 


P8W^m(l'); 


P8V'-m(l') 




1 


8mm 


8221' 


(r8,8j)(m,2;e)(e,e') 





1 

2 


P8W*m(l'); 


P8^b^'m{l') 










(r8,8|)(m,2,-)(e,e') 





1 
? 


P8^'m2^m(l'); 


P?,^'b'^m{l') 


222 


4mm 


2 


4*22* 


(r8,4|)(m,e)(e,2,)(e,2y-) 


A 


f8'**m«m(222); 


P%'^' bam{222) 




4 


2mm 


4221' 


(r8,4,-)(m,e')(e,2;E)(e,2y-) 


A 


PS*m'„m{222y, 


PS%m{222) 
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Table 6 

Restrictions on the form of S(k) for any wave vector k in the magnetic lattice when Fe = 2, 2', or 1'. In each case the form of S(k) depends on the parity of J2 "i 
where k = ^ «, b'''. Colors refer to the points in Figure 3. 





"i even 


Hi odd 


Te 


Red and Black 


Green and Blue 


2 


(0,0,5,) 




2' 


(Sv,5,,0) 


(0.0, 5j) 


1' 


(0,0,0) 





Table 7 

Restrictions on the form of S(k) for any wave vector k in the magnetic lattice when Ye = 222 or 22'2'. In each case the form of S(k) depends on the parities of 
n\ + nj, and n2 + nn, where k = ^ Kjb''' . Colors refer to the points in Figure 3. 





ni + ns even 


ni + Hi odd 


ni + n-i odd 


ni + Hi even 




«2 + "4 even 


112 + "4 odd 


n2 + n4 even 


n2 + n4 odd 




Red 


Black 


Green 


Blue 


222 


(0,0,0) 


(0,0,5,) 


(5., 0,0) 


(0,5„0) 


ll'T! 


(0,0,5.) 


(0,0,0) 


(0,S„0) 





Table 8 

Additional restrictions on the form of S(k) for special wave vectors that are invariant under mirror reflections when G = 8mm. Note that in most cases the group of 
spin-space rotations F is abelian (except when F contains a 4-fold or an 8-fold rotation, with optional primes) in which case ' = ™d <5'/i5'^' = 5^, 

and the selection rules take a much simpler form. Vectors k, along mirrors m, with «, odd, and vectors k- along mirrors m- with «, + U odd, are represented as open 
circles in Figure 3. 



Spin space- 
group type 


k,- = n,-b(') +/i(b''-') -t-b('+')) 
along m, 


k 


; = «,(b<'' + b('+')) + /;(b('-i) + b(-+2)) 
along m- 




5'-V<5'"'S(k,-) =S(k,-) 


s'ns'-'s{ki) = s(k;) 




1 


5'-'^i5'-'S{ki) = S(k,) if n; even 
S'-^^iS'-'S{ki) = -S(k,) if n; odd 


5'M<5>-'s(k;) = s(k;) 




5'-V'5'"'S(ki) = S(ki) 


1 


<5'Aj5'-'S(k,') = S(k,') if Hi + /,■ even 
, <5'M5'-'S(k,') = -S(k,') if ra,- + /; odd 




1 


<5'- 1 -'S(ki) = S(ki) if m even 
, 5'- V5'"'S(ki) = -S(k,) if n,- odd 


1 


<5'/x<5i-'S(k;) = S(k;) if n; + /; even 
, S'uS^-'Sm) = -S(k;) if n; + /; odd 
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Choice of L from a known Bravais class 

i ^ 
Choice of G 

i 

Choice of G, 



i 





Calculation of possible F and 






i 






Pairing of cosets 






i 






Spin < 


,eomeli'ic cryslal c 


lasses 






i 




Choice of relative orientation of L and G 




i 






Spin arithmetic crystal classes 





i 



Choice of gauge 

i 

Solution of group compatibility conditions 

i 

Spin space-group types 



Figure 1 

Flowchart describing the steps required for tlie enumeration of spin point groups and spin space-groups. Double boxes indicate the classification into equivalence 
classes as described in section 3. 




Figure 2 

Generating vectors and mirror lines for two-dimensional octagonal lattices. The solid arrows are the star of generating vectors and their negatives ±bW 
The dashed lines show the two types of mirrors in the 8mm point group, as described in the text. 



.±bW. 
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Figure 3 

(Color) A subset of the wave vectors of the 2-dimensional octagonal lattice, exhibiting all possible selection rules. The lattice generating vectors b''' and their 
negatives, as well as the origin, are denoted by solid black circles. The rest of the points shown are of the form k = nib''' + n2b'^' + nsb'^' + n4b''*', with indices 
running from —6 to 6. Colors encode the parities of the indices of k at each point as follows. Red: iii + and n2 + rn both even; Black: n[ + and n2 + «4 both 
odd; Blue: n[ + rij even and n2 + 114 odd; Green: n[ + 113 odd and fi? + 'U even. These color codes shold be used together with Tables 6 and 7 to determine the 
selection rules at each wave vector that are due to the lattice spin group Fj. Vectors k, = fi,b''' + /,(b''~'' + b''+'') invariant under mirrors m, with n, odd, and 
vectors k' = n, (b''' + b''^'') + Z, (b''~'' + b'''*'^') invariant under mirrors m'- with «, + /, odd, are represented as open circles. These points should be used together 
with Table 8 in determining the additional selection rules for wave vectors along mirror lines, when the point group is 8mm. 



